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After proving and using a mathematical identity, we first deduce an expansion formula of operator 
binomials power. Then, starting from our idea of combining the Feynman path integral spirit and 
the Dyson series kernel, we find an explicit and general form of time evolution operator that is 
a c-number function and a power series of perturbation including all order approximations in the 
unperturbed Hamiltonian representation. Based on it, we obtain an exact solution of the Schrodinger 
equation in general quantum systems independent of time. In special, we write down the concrete 
form of the exact solution of the Schrodinger equation when the solvable part (unperturbed part) of 
this Hamiltonian is simply taken as the kinetic energy term. Comparison of our exact solution with 
the existed perturbation theory makes some features and significance of our exact solution clear. 
The conclusions expressly indicate that our exact solution is obviously consistent with the usual 
time-independent perturbation theory at any order approximation, it explicitly calculates out the 
expanding coefficients of the unperturbed state in the non-perturbation method, and it fully solves 
the recurrence equation of the expansion coefficients of final state in the unperturbed Hamiltonian 
representation from a view of time-dependent perturbation theory. At the same time, the exact 
solution of the von Neumann equation is also given. Our results can be thought of as theoretical 
developments of quantum dynamics, and are helpful for understanding the dynamical behavior and 
related subjects of general quantum systems in both theory and application. Our exact solution, 
together with its sequence studies on perturbation theory [An Min Wang, quant-ph/0611217] and 
open system dynamics [An Min Wang, quant-ph/0601051] can be used to establish the foundation 
of theoretical formulism of quantum mechanics in general quantum systems. Further applications 
of our exact solution to quantum theory can be expected. 

PACS numbers: 03.65.-w, 03.65.Ca 



I. INTRODUCTION 

One of the most important tasks of physics is to obtain the time evolution law and form of a physical system, that 
is, so-called dynamical equation and its solution. This directly relates to physics foundation. In a quantum system, 
the dynamical equation is the Schrodinger equation [H, 0] for a pure state or the von Neumann equation for a 

mixed state p{t) [Ij, that is 

-i^|M/(i)) = i/|*(t)), (1) 



p{t) = -nH,p{t)]. (2) 



where H is the Hamiltonian of a quantum system. However, only a few idealized quantum systems are exactly solvable 
at present by using the existed theory and methods. Therefore, it is extremely interesting and essential important 
to find the exact solutions to the Schrodinger equation and the von Neumann equation for Hamiltonians of even 
moderate complexity. 

If the Hamiltonian is assumed independent of time, an arbitrary initial state at i = is denoted by |^(0)) for a 
pure state or p{0) for a mixed state, then the final state at a given time t will become 

|vl/(i))=e-'«*|*(0)), (3) 
p{t) ^ e-'^*p(0)e'^*, (4) 

where e~''^* (set /i — 1 for simplicity) is called as the time evolution operator. The above equations indicate that 
that the final state can be rewritten as the form that the time evolution operator acts on an initial state. In other 
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words, they are, respectively, a formal solution of the Schrodinger equation for a pure state and one of the von 
Neumann equation for a mixed state. However, the formal solution is not the practical solution that can be applied to 
concrete problems and actual calculations. Ones are interested in the practical solution because it can really produce 
the physical results and conclusions in the quantum theory. Hence, to solve the Schrodinger equation and the von 
Neumann equation refers to find the practical solution in most cases. It is clear that this task is equivalently to seek 
for the explicit expression of the time evolution operator. 

Historically, the most famous example to study the expression of the time evolution operator is Feynman path 
integral formulism In fact, Dyson series is also an important example 5]. From our point of view, only if the 
expression of the time evolution operator is a c-number function, can we clearly express the practical solution of the 
final state using ^ and ([U, respectively, for the pure state and the mixed state. Moreover, we realize that for a 
general quantum system without the exact solution in the usual theory, the expression of the final state in c-number 
function form can be generally expressed by an infinite series. Hence, only if the expressing series of the final state 
is a power series of perturbation, can we inherit the advantages of perturbation theory. As is well-known, to find the 
c-number function form of the time evolution operator is a successful linchpin of Feynman path integral formulism, 
and to expand the time evolution operator as a power series of perturbation is a powerful headstream of Dyson series 
(in the interaction picture) . Consequently, our physical idea in this paper comes from the combination of Feynman 
path integral spirit and Dyson series kernel. Starting from this combined idea, we would like to find an explicit and 
general expression of the time evolution operator that it not only is a c-number function but also is a power series 
expansion of perturbation including all order approximations in terms of our own methods. 

To express the time evolution operator as a c-number function needs a representation of Hilbert space, and further 
to expand it as a power series of perturbation needs to separate its unperturbed part and perturbing part. Actually, 
the above problems can be studied by using the similar methods to the usual perturbation theory [l|, y] . 

In the usual perturbation theory, the key idea to research the time evolution of general quantum systems is to split 
the system Hamiltonian into two parts, that is 

H = Ho + Hi, (5) 

where the eigenvalue problem of so-called unperturbed Hamiltonian Hq is solvable, and so-called perturbing Hamil- 
tonian Hi is the rest part of the Hamiltonian. In other words, this splitting is chosen in such a manner that the 
solutions of Hq are known as 

Fo|*^)=i?7l*^): (6) 

where 1$''') is the eigenvector of Hq and E-y is the corresponding eigenvalue. Whole \^^), in which 7 takes over all 
possible values, form a representation of the unperturbed Hamiltonian. It must be pointed out that in the perturbation 
theory, the above Hamiltonian split principle with the best solvability is not unique in more general cases. It is possible 
to need to consider such an additional condition that the selected split can remove degeneracies as complete as possible 
in order to suit to apply the perturbation theory, or to let our improved scheme of perturbation theory Q work well. 
If the remained degeneracies are allowed, it requires that the off-diagonal element of the perturbing Hamiltonian 
matrix between arbitrary two degenerate levels are always vanishing. As an example, it has been discussed in our 
sequence study Q. Of course, if the cut-off approximation of perturbation is necessary in our improved scheme of 
perturbation theory Q , the off-diagonal elements of Hi should be small enough compared with the diagonal elements 
oi H = Hq + Hi in the Hq representation. 

As soon as we find the explicit and general expression of the time evolution operator that is both a c-number 
function and a power series of perturbation including all order approximations, we definitely can obtain the exact 
solutions of the Schrodinger equation and von Neumann equation, and so we can solve a lot of problems in quantum 
theory, for example, entanglement dynamics and open system dynamics. Of course, our solutions are called "exact 
ones" in the sense including all order approximations of perturbation. However, our final purpose is not limit to these. 
We, starting from our exact solutions, try to finally establish the foundations of theoretical formulism of quantum 
mechanics in general quantum systems. Our serial studies on perturbation theory Q and open system dynamics Q 
indicate that we have partially arrived at our purpose. 

In this paper, the key and central task is, in the unperturbed Hamiltonian representation, to find an explicit 
expression of the time evolution operator e~'^' that is a c-number function and a power series of perturbation 
including all order approximations according to our idea combining the Feynman path integral spirit and Dyson 
series kernel. Our method is first to derive out an expansion formula of operator binomials power and then apply 
it to the Taylor's expansion of time evolution operator e~'^* = e~'(^"+^i)*. Moreover, after proving and using our 
identity, we derive out the explicit and general form of representation matrix of e~'^* in the representation of the 
unperturbed Hamiltonian Hq. Consequently, we obtain the exact solutions of the Schrodinger equation and the von 
Neumann equation in general quantum systems independent of time, in particular, a concrete form of the Schrodinger 
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equation solution when the solvable part of Hamiltonian is simply taken as the kinetic energy term in order to 
account for the generality and universality of our exact solutions. Furthermore, by comparing our exact solution of 
the Schrodinger equation with the usual perturbation theory, we reveal their relations and show what is more and 
what is different. The conclusions clearly indicate that our exact solution is obviously consistent with the usual time- 
independent perturbation theory at any order approximation, but also in our exact solution we explicitly calculate out 
the expanding coefficients of unperturbed state in Lippmann-Schwinger equation |8j for the non-pcrturbative method, 
and/or fully solves the recurrence equation of the expansion coefficients of final state in the unperturbed Hamiltonian 
representation from a view of time-dependent perturbation theory. Based on the above all of reasons and results, our 
exact solutions can be thought of as the theoretical developments of quantum dynamics in general quatum systems, 
and they are helpful for understanding the dynamical behavior and related subjects of quantum systems in both 
theory and application. From our point of view, our exact solutions and their sequence studies on perturbation theory 
and open system dynamics can be used to establish the foundation of quantum mechanics of general quantum 
systems. Specially, we think that the features and advantages of our exact solutions can not be fully revealed only 
by the improved scheme of perturbation theory @ and open system dynamics We would like to study the more 
applications to the formulation of quantum mechanics in general quantum systems in the near future. 

This paper is organized as the following: in this section, we give an introduction; in Sec. |TT]we first propose the 
expansion formula of power of operator binomials; in Sec. Illll we derive out an explicit and general expression of time 
evolution operator after proving and using our identity; in Sec. IIVI we obtain the exact solutions of the Schrodinger 
equation and the von Neumann equation, and then present a concrete example when the solvable part (unperturbed 
part) of Hamiltonian is taken as the kinetic energy term; in Sec. |V] by comparison of our exact solution of the 
Schrodinger equation with the usual perturbation theory, we prove their consistency and their relations, and explain 
what is more and what is different.; in Sec. IVII we summarize our conclusions and give some discussions. Finally, we 
write an appendix where the proof of our identity is presented. 



II. EXPANSION FORMULA OF OPERATOR BINOMIALS POWER 



In order to obtain the explicitly exact solutions of the Schrodinger equation and von Neumann equation in general 
quantum systems independent of time, we need to deduce the expression of the time evolution operator. According to 
our physical idea, we should first separate the time evolution operator into two parts, respectively, with perturbation 
and without perturbation. After to split the Hamiltonian, this problem changes to the derivation of expansion formula 
of operator binomials power. Without loss of generality, we are always able to write the power of operator binomials 
as two parts 

(^-KB)" = + (7) 

where A and B are two operators and do not commute with each other in general. If B is taken as a perturbation, 
the above decomposition just splits the operator binomials power into two parts, respectively, with perturbation and 
without perturbation. Now, the key matter is how to obtain the form of /"(^, B). 

It is clear that /"(A, B) is a polynomial including at least first power of B and at most nth power of B in every 

term. Thus, a general term with Ith power of B has the form ^01=1 A'^'B^ fe; ^ From the symmetry of 

power of binomials we conclude that every ki take the values from to (n — /), but it must keep n — l — X]i=i ^ 0. 
So we have 



n n — l /I \ 

r\A,B) = Y^ rQA'=.Br4"-'-^'=i'=' (s) 

1=1 fci,--- ,fci=0 Vi=l / 

-EE ^l[^'^B]A--^-^'--^'^^9in-l-Y.k^■ (9) 

;=1 fci, -- ,fc;=0 \i=l / V i=l ) 

where B{x) is a step function, that is, d(x) = 1 if a; > 0, and d(x) = if a; < 0. Obviously based on above definition, 
we easily verify 

f{A,B)=B, f{A,B)=AB + B{A + B). (10) 



They imply that the expression ([5]) is correct for n = 1,2. 
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Now we use the mathematical induction to prove the expression ^ of f'^{A,B), that is, let us assume that it is 
valid for a given n, and then prove that it is also valid for n + 1. Denoting B) with the form of expression 

([9]) where n is replaced by n + 1, and we extract the part of / = 1 in its finite summation for I 

We extract the terms fci = in the first and second summations, and again replace fci by fci — 1 in the summation 
(the summations for ki from 1 to (n + 1) — / change as one from to n — Z), the result is 

n n+l {n+l)-l / ' \ / ' \ 

= BA"" +Y A^'BA'^-^' +bY^ myl'^'B L4("+i)-'-^'=2'^'6' (n + l)-/-^fc, 

fci=l i=2 fc2, - ,fei=0 \i=2 / \ z=2 / 

„+l(n+l)-i (n+l)-i / ' \ / ' \ 

+ E E E [\{A'^B]A^-^^^-^-^'^^^^9[{n + l)-l-YkX (12) 



furthermore 



1=2 fci = l fc2,---,fci=0 \i=l / \ 1=1 / 

n-1 n+l (n+l)-i / I \ 

r''+^{A,B) ^ BA'' +AYa'''BA"-^-^' +bY E n^'''^ 

fci=0 i=2 fc2r-- ,fcl=0 \i=2 / 

X A("+i)-'-5:t2 (^(r^ + 1) - ; - ^ fc.^ 

n+l n-i (n+l)-i / ' \ / ' \ 

+^EE E n^''^P""'"^-^''^h-'-EfcO- (13) 

i=2 fci=0 A:2,- - ,fcl=0 \i=l / V i=l / 

Considering the third term in above expression, we change the dummy index {A;2, ^3, • • ■ , ki\ into {fci, fc2, • ■ ■ , h^i}, 
rewrite (n + l) — / asn — (/ — 1), and finally replace Z by Z — 1 in the summation (the summation for I from 2 to n + 1 
changes as one from 1 to n), we obtain 

n+l (n+l)-i f I \ f ^ \ 

^E E n^''^p("+i)-'-^'-^'=»0 (n+l)-/-^fcj 

i=2 fc2r - .fc!=0 \i=2 / V 1=2 / 

n+l n-(i-l) n-\ \ / \ 

= ^E E n^'"^P"-('-i)-^'=5'='(? n-(?-l)-^fc, 

i=2 fci, -- ,fci_i=0 \i=l / \ i=l / 



n n 



-I 



-^E E n^^-'^P"-"^'-^-'^ --'-E^O 

i = l fci, -- ,A;,=0 \i=l / \ i=l / 

^BriA,B), (14) 

where we have used the expression (HI) of /"(A, B). Because of the step function 9 ( n — I — X]i=i i'^ the fourth 

term of the expression (jl3p . the upper bound of summation for Z is abated to n, and the upper bound of summation 
for A;2, fcs, • ■ • , fc; is abated to n — I. Then, merging the fourth term and the second term in Ea. p3|) gives 



n-l n+l n-i (n+l)-i / ' \ / ' \ 

aYa'^^ba^-^-^^+aYY. E n^''^M""'"^'^'''M"-^-E^O 

fci^O 1^2 ki^O k2,--- ) \ %^\ I 

n— 1 n n—l / I \ / ^\ 

^aya'^'ba^'^-^^+ay y n^'"^p"~'"^'^^'"^ "-^-E'^0 

fci=0 (=2 fci,fc2,--- ,fci=0 \i=l / V i=l ) 

n n~l / ' \ / ' \ 

= ^E E n^'"^P""'"^'^^''M"-'-E'^0 

(=1 fei,fc2,- - ,fci=0 \i=l / V i=l / 

= Ar(A,S), (15) 
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where we have used again the expression ([9]) of f"{A, B). 

Substituting ([H]) and (|15p into (fT^]) . immediately leads to the following result 



Note that 



we have the relation 



jr«+i = BA" + {A + B)f'\A, B) 



{A + BY+^ ^ {A + B){A + BY 



B) = BA'^ + {A + B)r{A, B) 



Therefore, in terms of Egs. pB]) and (fT5)) we have finished our proof that f'^{A^B) has the expression 
any n. 



(16) 



(17) 



(18) 
or © for 



III. EXPRESSION OF THE TIME EVOLUTION OPERATOR 

Now we investigate the expression of the time evolution operator by means of above expansion formula of power of 
operator binomials, that is, we write 



iHt 



n=0 ' ri=0 



(19) 



In above equation, inserting the complete relation !$''')($'''[ — 1 before every Hq\ and using the eigen equation 
of Hq ([S]), it is easy to see that 

I 



f 



n n — l V I 

i=l 7i,--- ,7i + l fc,.---.fc,=0 U=l 



i=l 7i,--- ,7i + l fci,---,fci=0 



7i 



-^^7! + ! 



.2=1 



(j)7l|)^(j)7! + l| 



|$7l)($7! + l|^ 



1=1 7i,--- ,7i+l U=l 

where H^^''^''^^ = ($^' £^[7, ^] is a vector with / + 1 components denoted by 

I] = {^^71 , £^72 , • ■ • , £"7, , £^7,+i } 

and we introduce the definition of Cl''{Eh, I]) {I > 1) as the following 



cr{E[j,i])^ 

5^^^^ ki+l<; 



.1=1 



^71 + 1 



(20) 



(21) 



(22) 



The expression (j20p and the expression of the time evolution operator depend on it has incarnated our physical idea, 
that is, it is both a c-number function and a power series of perturbation. Our further task is to find the explicit 
expression of C"(£'[7, 1]). 

Before to find C"{E['-f,l]), let us mention the cases when energy level degenerations happen, that is, some same 
values appear in the eigen spectrum (£ = {-B7I7 = 1, 2, ■ • • }, the eigen spectrum is re-denoted by € = {Eja^^\j = 
2. • ■ • . a.,, = 1. • • ■ , m-y} and the numbers of original 7 is different from the numbers of new 7. Here, m-y is called 



1,2.--- = 1 

ijij — 1, the energy level E^a-y is not degenerate 



, il, • ■ ■ , Lt^ ±, • • ■ , lll-yf dllU. LllC ilLllllUCia Ul K^Ll^lliai -f 10 UIliCiCllL iiUlll LilC llUillUCiD Ul new - J. AACiC, / / t-y lO L-CLllCU 

degeneracy for a given 7, which means the number of the same eigenvalues (m^-fold degeneracy). In special, when 
' " ' ' It is clear that the above expression (PP]) has only a little change: 



n 

nHo,Hi) = E E E Q"(^[7:^]) 

Z=l 7l,--- >7i + l Oti ^o-ji^i 



.1=1 



|$7ia7i)($'''' + i°"'i + l 



(23) 
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which has no obvious influence on our following derivation and proof. For simplicity, we do not consider the degenerate 
cases in this section since no new idea and skill will be needed here. 

In order to derive out an explicit and useful expression of C"(i?[7, Z]), we first change the dummy index ki —> 
n — I — X]i=i ill Eq.®. Note that at the same time in spite of 9{n — I — X]i=i ^i) changes as 0{ki), but a hiding 
factor 9{ki) becomes 9{n — / — ^O- Thus, we can rewrite 



n n 



-I /l-l 



(24) 



1=1 ki.--- ,ki=0 \i=l 



In fact, this new expression of f^{A,B) is a result of the symmetry of power of binomials for its every factor. 
Similarly, in terms of the above method to obtain the definition of C"{E[j, I]) (|22p . we have its new definition (where 
ki is replaced by k) 



n — l I n—l 



cr(i?[7,/]) = E E 



fc=0 I fci,--- ,fci_i=0 U=l 



7-1 



i-i 



E\r''-^-^''-'0[{n-k)-l-J2kA\E'^^^^. (25) 



Because ^(n — k) — I — X]i=i ^* j abates the upper bound of summation for ki{i = 1, 2, — 1) from (n — /) to 
(n — fc — 1) — (Z — 1), we obtain the recurrence equation 



f (n-fc-l)-(i-l) 



cr{E[jj]) = j2{ E 

k=0 [ fci,--- ,fc,_i=0 



1-1 



,{n-k-l)-j:'rlki~{l-l) 



i=l 

xe (^{n-k-i)-{i- 1) - E^») } <+i 

n — l 

Ecr^i(i?[7,z-i]X,,. 



fc=0 

In particular, when I = 1, from the definition of Cl'{E{'-f, 1)) it follows that 

n—l n—l 



fci=0 fei=0 

By means of the summation formula of a geometric series, we find that 



cr(i?(7,i)) = E KiK-'-''-^--' E ( 



Based on the recurrence equation (j26[) . we have 



EZ 



77' n 

^72 



E — E E — E 

^71 -'^72 ^71 -^^72 



C^Hi?(7,2)) = 



n-2 / -^ri-k-1 



Tpn—k— 1 

^72 



fe=0 
2 



J^'-y-^ -Ej-yey -Ej-yey 

-1 z^n-l n-2 ^ X fc 



Et 



\p (-1)' -E"; X - / £^ 



k=0 



{^iy-\El - E^^El"') 
~{ (E-n - E^^){E^^ - Ej^) 



e: 



71 



e: 



72 



(^71 ^ Ey^){E^^ - Ey^) {E^^ - Ej^){E^2 - E^^) 

^ fSa 

(^71 ~ Ej^){Ej2 ~ -£'73) 



(26) 



(27) 



(28) 



(29) 
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In the above calculations, the last step is important. In fact, in order to obtain the concrete expression of C" {I and 
n are both positive integers), we need our identity 

'+1 f (If < if < I) 

It is proved in Appendix A in detail. The denominators di{E['y, I]) in above identity are defined by 

di{E[-f,l]) - n(^7i~^7.+i). (31) 

2=1 

d^{Eb, ;]) - n (^7. - ^7,) n (^7. -E-,J, (32) 

j=l k=i+l 

I 

di+,{E[^,l]) - n(^7. "^7,+i). (33) 
where 2 < i < I. Then, using the recurrence equation (pSj) and our identity (PU)) . we obtain 

Here, the mathematical induction shows its power again. If the expression (IM)) is correct for a given n, for example 
n = 1,2, then for n + 1 from the recurrence equation (p6)) it follows that 



^^dmi,i-i]) to V ^7 

t d.(i?[7,/-l])(i5^.-i?^,+J ■ ^ ' 

Because that di{E[y, I — 1])(E^. — E^^^^) = di{E[j, I]) {i < I + 1) and our identity pO|) . it is easy to see 

MEb,i]) ^ ' di+,{E[^,i]y ^ ' 



cr-(i^[7,^]) = E(-ir^^^^- (37) 



Substitute it into Eq. (|35|) yields 



2 = 1 



That is, we have proved that the expression of CJ^ {E[y , I]) is valid for any n. 

It must be emphasized that since the summation is over all of values of the energy eigenvalues, there are ap- 
parent divergences in the expression of Cp{E['y,l]). Here, "apparent" refers to an unture thing, that is, the ap- 
parent divergences are not real singularities and they can be eliminated by mathematical and/or physical meth- 
ods. Moreover, the same problem appears if with degeneracy. Therefore, we need to understand above expres- 
sions in the sense of limitations. For instance, for C"(i?[7, 1]) — {El!j_^ — i?"^) / [E^^ — E^^), we have the expression 
limE^^^E.,-^ C^'iEl-f,!]) = {5jij2 + Q{\j2 - 1i\)Se^-^e^J nE:^-^ , where the step function e(x) = 1, if a; > 0, and 
&{x) = 0, if X < 0. In other words, the apparent divergences here are not real divergences and can be eliminated by 
finding the correct limitations. There is no theoretical problem when we formally keep the apparent divergences in 
the above expression. Actually, in our recent work [6|, we successfully eliminate all the apparent divergences. 
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Because of our identity ((30)) . the summation to I in Eq. (|20)) can be extended to oo. Thus, the expression of the 
time evolution operator is changed to a summation according to the order (or power) of the perturbing Hamiltonian 
Hi as foUows 



) e 



1 = 1 ll,--- ,71 + 1 



d,{E[j,l]) 



ni:r7j7j + ir r 
111 "7l7"7l + n 



When there is degeneracy, it becomes 



($^|e-' 



Ht 



$7 



G S-y '•y ' 



EE E 

1=1 7l>-- - >7i + l jO^,^! 



'l+l 

E(-ir^ 

.i=l 



d,(£;[7,;]) 



j=i 



(=0 



$7 



that is. 



where 



Z=0 Z=0 7,7' 



E 



7l,--- ,11 + 1 \-i=l 



'l+l 



iE^.t 



J-J^7,7. + l 

3 = 1 



'^7l7'^7i + l7'- 



(38) 



(39) 
(40) 

(41) 

(42) 
(43) 



It is clear that this propagator has the closed time evolution factors. In special, it is a c-number function just 
like the Feynman path integral Q, and it is a power series of the perturbing Hamiltonian just like Dyson series in 
interaction picture [5]. Therefore it has, at the same time, the advantages that exist in the Feynman path integral 
famulism and Feynman diagram (Dyson) expansion formula . 



IV. SOLUTION OF THE SCHRODINGER EQUATION AND THE VON NEUMANN EQUATION 

Substituting the expression of the time evolution operator (j38p into Eq. ([3]) , we obtain immediately our explicit form 
of time evolution of an arbitrary initial state in a general quantum system 



\^{t)) ^^,(t)|vi/(o)) = ^^A7^'(i)[($V|vi/(o))]|$^), 

1=0 7,7' 



(=0 



(44) 



In form, this solution is also suitable for the degenerate cases, but more apparent divergences will appear. It must be 
emphasized that these apparent divergences can be easily eliminated by the limitation calculations. To save space, we 
delay to discuss the disposal of degenerate cases and the elimination of apparent divergences to our sequence study 

Similarly, the exact solution of the von Neumann equation for mixed states can be obtained by using our general 
and explicit expression of the time evolution operator (|4ip . From the formal solution of the von Neumann equation 



p{t) = e-'^V(0)e'^* 



it immediately follows that 



Pit) 



£ Ak{t)piO)Ai{-t) £ Ak{t)piO)AUt) 

oc 

E EE<w/''(oMf'''(-oi'!>'^)(<i'^' 



k,l=0 7,7' 13,0' 



(45) 

(46) 
(47) 
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Since wc have obtained the exact expression of the time evolution operator, it is direct to extend our conclusions 
from the pure state to the mixed state. For simplicity, we only study the cases of pure states in the following, that is, 
our exact solution of the Schrdinger equation. 

Note that the linearity of the evolution operator and the completeness of eigenvector set of Hq, we can simply set the 
initial state as an eigenvector of Hq. Up to now, everything is exact and no any approximation enters. Therefore, it is 
an exact solution despite that its form is an infinity series. In other words, it exactly includes all order approximations 
of Hi . Only when Hi is taken as a perturbation, can it be cut off based on the needed precision. From a view of 
formalized theory, it is explicit and general, but is not compact. Moreover, if the convergence is guaranteed, it is 
strict since its general term is known. Usually, to a practical purpose, if only including the finite (often low) order 
approximation of Hi , above expression is cut off to the finite terms and becomes a perturbed solution. 

Actually, since our solution includes the contributions from all order approximations of the rest part of Hamiltonian 
except for the solvable part Hq, or all order approximations of the perturbative part Hi of Hamiltonian, it is not 
very important whether Hi is (relatively) large or small when compared with Hq. In principle, for a general quantum 
system with the normal form Hamiltonian, we can always write down 



H = 



2m 



V. 



(48) 



If we take the solvable kinetic energy term as Hq 



2m 



and the potential energy part as Hi = V, our method is 



applicable to such quantum system. It should be pointed out that we study this concrete example in order to account 
for the generality and universality of our exact solution. 



It is clear that 



And denoting 



2m 



{x\k) 

MO) 

*(a;,t) 



1 



Akx 



mm, 



(49) 



(50) 

(51) 
(52) 



we obtain the final state as 



k /=! fe,fe' fei ,fei4_i 



l+l 



{-ly 



^^ di{E[k,l]) L3/2 



i(fe-a!— -Efcjt) 



(53) 



i(fcx-Efct) .^.^ 

E^^(o) +EE E v-fe'W 
fe (=1 fe,fe' fei ,fei-i-i 



nv(fe,-fe,+i) 



E 



J_ i(fe-a!-Bfc.t) 



f^di{E[k,l]) £3/2 

In the second equal mark, we have used the fact that V = V{x) usually, thus 

= {kj\V{x)\kj+i) 

/oo 
d^Xjd^Xj+i {kj\xj){xj\V{x)\xj+i){xj+i |fej+i) 
-OO 

= -3 / d'xjd'xj+ie-"^r^^V{xj+i)SHxj - x,+i)e''=^+--+^^- 
= -3 / d3xy(a;)e-'('=^-'=^+i>'= 

J —OO 

= V(fcj - kj+i). 



(54) 



(55) 
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It is clear that V(fc) is the Fourier transformation of V{x). Eq. (|53)l or ((54| is a concrete form of our solution (|44|) 
of Schrodinger equation when the solvable part of Hamiltonian is taken as the kinetic energy term. Therefore, in 
principle, if the Fourier transformation of V{x) can be found, the evolution of the arbitrarily initial state with time 

can be obtained. It must be emphasized that it is often that the solvable part of Hamiltonian will be not only T 

for the practical purposes. 

Furthermore, we can derive out the propagator 



2m 



X') 



{x\k){k\e~'"'\k')(k'\x') 

k,k' 

oc 



l—l h,h' k\. 



■J~J ykjkj + i 



Sk^k^ki+^k' 



l + l 



E 



(-1)^ 



^^d,{E[k,l]) L3/2 



g-i-Efc.t 



ik-x — \k'-x' 



(56) 



It is different from the known propagator in the expressive form, but they should be equivalent in theory because we 
still use the usual quantum mechanics principle. However, what is more and what is new in our exact solution? In 
the following, and in our serial studies [1, 0| , we will see them. 



V. COMPARISON WITH THE USUAL PERTURBATION THEORY 



In this section, we will compare the usual perturbation theory with our solution, reveal their consistency and 
relation, and point out what is more in our solution and what is different between them. Furthermore, we expect to 
reveal the features, significance and possible applications of our solution in theory. We will respectively investigate 
and discuss the cases comparing with the time- independent, time-dependent perturbation theories and the non- 
perturbative solution method. 



A. Comparison with the time-independent perturbation theory 

The usual time-independent (stationary) perturbation theory is mainly to study the stationary wave equation in 
order to obtain the perturbative energies and perturbative states. However, our solution focus on the development of 
quantum states, which is a solution of the Schrodinger dynamical equation. Although their main purposes are different 
at their start points, but they are consistent. Moreover, our exact solution is also able to apply to the calculation of 
transition probability, perturbed energy and perturbed state, which will be seen clearly in Ref. Moreover, our 
exact solution has been successfully applied to open system dynamics Q- 

In order to verify that the time-independent perturbation theory is consistent with our solution at any order 
approximation, we suppose that in the initial state, the system is in the eigen state of the total Hamiltonian, that is 



(57) 



It is clear that 



l^EAt)) = e-^-*|vl/^,(0))=5] e-^-*V+^A7^'(0 [(<i>^>i^,(0))] 1$^) 

7,7' I ' = 1 J 

Calculating the Kth time derivative of this equation and then setting t — 0, we obtain 



(58) 



7,7' 



1=1 



(59) 
(60) 



7,7' 



11 



where we have used the fact that 



(61) 



It is easy to obtain that 



E 



71 ,7! + l 



'^7l7'^7! + l7' ' 



Br\K) = ^ Cf(i?[7,Z])0(i^-O 

7l,--- ,11 + 1 

where we have used our identity ((30)) . This means that \i I > K 



]-J^7.7, + l 

i=i 



'^7l7'^7l + l7' 



In special 



dt^ 



AAr it) 



{^i^Br {K>l) = Q. 



At 



= -iBr (1) = -inr Sn 



t=0 



(62) 



(63) 
(64) 



(65) 



(66) 



d^Aj^ jt) 



dt2 



(67) 



In other words, the summation to / in the right side of Eq. (|59)) is cut off to K. Therefore, by left multiplying ($^| to 
Eg. ([55)1 we obtain 



K 



Er'^a^ = E^a^ + E E 

7' 1=1 

When K = 1, the above equation backs to the start point of time- independent perturbation theory 

7' 



(68) 



(69) 



It implies that it is consistent with our solution. In order to certify that our solution is compatible with the time- 
independent perturbation theory, we have to consider the cases when K > 2. By using E^ — Ej, i^B^ ^a^^^ 
substituting Eq. (|68p twice, then moving the second term to the right side, we have 



K 



K-l 



E^-^Era^ = E^^a^ + E E {K)a^' ^ E E ^^'(^ " 1) E (^7"^7'7" + ^7'^ j «• 



7' 1=1 

K-l 



7' 1=1 



7" 



<«7 + E E [Br\K) EyB^'iK - 1) 



7' 1=1 
K 



K-l 



E [li^i'"''' '577iVi7'«7' -E E Br'iK-i)j2Hr"c 

7' \J = 1 / 7' '=1 7" 



(70) 
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Since when K > 2 

Cf (i?[7,l]) - E,,C^-\E[^,1]) = E^-\ (71) 

and again / > 2 

C^{E[^,l]) - E^,^,C^-\E[^,l]) = Cl-\E[^,l- 1]), (72) 
which has been proved in appendix A. Obviously, when if = 2, Eg. ((70)) becomes 



E^Era^ = E^a^ + E^Y^ H^'^'ay. (73) 

7' 

It can back to Eg. ([55]) . Similarly, when K >3, so do they. In fact, from Egs. (|71l72p . Eg. ([70)1 becomes 

E^-^Era^ = E^ a-, + E^-^ ^ H^' ay , (74) 

7' 

where we have used the fact 

Y,Y.[Bi'''iK)-EyBr\K~l)] «7'+E (fi^i"^""^" I '^77.'5. 



'^77l"7( + l7'°'7' 



7' '=2 7' 
K-1 



E 5: i?7-'(if -1)Y: m'^"a,„ = 0. (75) 
7' /—I 7" 



Its proof is not difficult because we can derive out 

K-l 



Y^Y. [Br\K)-E,,Br'{K-i] 



7' 1=2 

EE E [Cf{Eh,l])-EyCr\E[^iM [i{Hr'^^\6,,,5. 

7' Z=2 7l.--- ,71 + 1 \> = ^ 

EE E C,^T\£^[7,^-l])(ni^?"^"M^77lW7'«7' 



771 "7i + l7 ""7 



= EE E c^r^(i^[7,^])(riffi™M^r'''<57 

7' 1=1 7l>--->7! + l = l / 

= E E^7''(^- W'''""^"- (76) 
7', 7' '=1 

The first equality has used the definitions of (K), the second equality has used Eq. ([72]) . the third equality sets 

I — 1 I and sums over 7/+2, the forth equality sets 7' 7", inserts (57, ^^^7' and uses the definition of B'^''' {K) 

again. Substituting it into the left side of Eg. (175)) and using B]^_^{K — 1) expression, we can finish the proof of 
Eg. (j75p . Actually, the above proof further verifies the correctness of our solution from the usual perturbation theory. 

In our point of view, only if the expression of the high order approximation has been obtained, can we consider 
its contribution in the time-independent perturbation theory since its method is to find a given order approximation. 
However, this task needs to solve a simultaneous equation system, it will be heavy for the enough high order ap- 
proximation. In Ref. ^6*1, we will give a method to find the improved forms of perturbed energy and perturbed state 
absorbing the partial contributions from the high order even all order approximations of perturbation. 
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B. Comparison with the non-perturbative method 



Now, let us see what is more in our solution than the nonperturbative method. Actually, since H is not explicitly 
time-dependent, its eigenvectors can be given formally by so-called Lippmann-Schwinger equations as follows [8|: 



l*2(±)) 



1 



-H,\^l{±)) 



1$'^) + - 
\^■^)+Gl{±)H^\^l{±)) 

bj^i — H ±iri 



(77) 
(78) 
(79) 
(80) 



where the complete Green's function G'^{±) = l/{E.y — H ±\rf) and the unperturbed Green's function G'q(±) = 
l/(-E^ — Hq ± ir;) satisfy the Dyson's equation 



G^(±) = Gl{±) + Gl{±)HiG'^ {±) 

OO 

= Y,{Gl{±)H,)' Gl{±). 



1=0 



Here, the subscript "S" means the stationary solution: 



(81) 
(82) 

(83) 



More strictly, we should use the "in" and "out" states to express it [9\. In historical literature, this solution is known 
as so-called non-perturbative one. It has played an important role in the formal scattering theory ^Sj. 
Back to our attempt, for a given initial state |5'(0)}, we have 



l*(0)) = ^(*2'(±)|*(0)>|<(±) 

7' 

Acting the time evolution operator on it, we immediately obtain 
l*(0) - E(*s'(±)l*(0))e-'^-'*|<(±) 

= E(*^'|E(^i^o'(T)) lvI'(0))e-'^v*^(G^'(±)i7j 



(84) 



7' 



E 

7,7' 



1=0 



l'=0 



$7 



V(*^Xo))+ E (<i>^'|(^i<(T)) |*(o))($^i(G;j'(±)ii, 



l,l'=0 



$7 



(85) 



comparing this result with our solution ([44]) . we can say that our solution has finished the calculations of explicit 
form of the expanding coefficients (matrix elements) in the second term of the above equation using our method and 
rearrange the resulting terms according with the power of elements of the perturbing Hamiltonian matrix H'l''^'*^ . 
This implies that our solution can have more physical contents and more explicit significance. Therefore, we think that 
our solution is a new development of the stationary perturbation theory. Seemingly, its physical results are consistent 
with the Feynman's diagram expansion of Dyson's series, but its general form is explicit, and convenient to calculate 
the time evolution of states with time. Specially, we can obtain the improved form of perturbed solution that absorbs 
partial contributions from the high order even all order approximations in our subsequent works [6], in which, it will 
be seen that our improved scheme of perturbation theory based on our exact solution has higher efficiency and better 
precision in the calculation. 



C. Comparison with the time-dependent perturbation theory 



The aim of our solution is similar to the time-dependent perturbation theory. But their methods are different. The 
usual time-dependent perturbation theory 2] considers a quantum state initially in the eigenvector of the Hamiltonian 
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where v^'^ (t) = 



Hq^ and a subsequent evolution of system caused by the apphcation of an exphcitly time-dependent potential V = 
Xv{t), where A is small enough to indicate y as a perturbation, that is, 

i^^\^it)}^{H, + Xv{tm{t)) (86) 

with the initial state 

l*(0)) - I*"). (87) 

It is often used to study the cases that the system returns to an eigenvector |<i>/) of the Hamiltonian Hq when the 
action of the perturbing potential becomes negligible. 

In order to compare it with our solution, let us first recall the time-dependent perturbation method. Note that the 
state |5'(t)) at time t can be expanded by the complete orthogonal system of the eigenvectors 1$''') of the Hamiltonian 
Hq, that is 

l*W)=E^7WI'i'^> (88) 

7 

with 

c^(t) = (<I>^|*(t)), (89) 

whereas the evolution equation ([55)1 will become 

d I 
ig^c-fit) = E^c^it) + X'^v'''^ (t)cY{t), (90) 

7' 

<^"''y Now setting 

c^{t) ^ b^{t)e-'^^\ (91) 



i|6,(t) = AEe'(^^-^.')V^'(t)6y(i). (92) 

7' 

Then, making a series expansion of b^{t) according to the power of A: 

6,(t) = 5]A'6«(t), (93) 

and setting equal the coefficients of A' on the both sides of the equation ([M]) . ones find: 

^§-/'Ht) = 0, (94) 

7' 

From the initial condition ([87)) . it follows that c-y(O) — 6^(0) — S^a — b'^\(i) and ^^■'(O) — if ^ > 1. By integration 
over the variable t this will yield 

= s^a, (96) 

&W(t) ^ _i rdiie'(^--^")*iu^"(ii), (97) 
Jo 

•^0 72 

= -J2 [ di2 / 'dtie'(-^^-^^2)t2ei(S-,2-£;")ti„772(^2)„72a(^^)^ (98) 



and inserting this into ([90]) give 



72 
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If we take v independent of time, the above method in the time-dependent perturbation theory is still valid. Thus, 

6W(t) = (e-'^°* - e-'^-*) e'^-*^;^", (99) 



72 



1 1 g— iB^jtgiiS^t 



1 

— Ea){Ey2 ~ Ea) 



72 



+ - _e-i^^=te'^-* 



(100) 



This means that 



2 

4'^*) - E E(-l)'"' rf^(J^^^^\]) ^^^^^'^77.'?72a, (101) 



W= E E(-1)" ^.,M^ ^.. e-'^^^*^,,,J,3„^;^^^-^;^-^^ (102) 



71,72,73 i=l ^ L/J 



Now, we use the mathematical induction to prove 



— e '^Ti ' — 



Obviously, we have seen that it is valid for I = 1,2. Suppose it is also valid for a given n, thus form Eqs. (|95p and 
h^Jj\0) = 0, ; > 1 it follows that 

ll;+^\t) = -i Tdr Ve'(^''-^-)^i;'^^6^")(T). (104) 
^0 ^ 

Substitute Eqs.jgi]) and (fT03)) we obtain 

Jo ^ 



/ " ^ ^' d^lE[j,n]) , "77l'^7.+ia 

" 7 7l,--- ,7,1+1 i=l ^ ^ ' -'^ \j=l 

^ '^±3 ('p-KiJ7,--ET„+2)* _ n / " \ 

%..5..g'-"" ^.l^h.ni)(^-^J., tn"^""J 

In the last equality, we have inserted '^7„+2/3i summed over 7 and integral over r. In terms of the definition of 

di{E[j, n], we know that di{E['y, n]){Ej. — E^^^^) = ^^(£'[7, n + 1]). Then based on our identity (|5D|). we have 

n+l 

S"(-IY-^ = -(-i)(«+2)-i 1 (106) 

' dmi.n+l]) ^ ^ d„+2(i?[7,n+l]) ^ ^ 
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Thus Ea. (|105p becomes 



"+2 g-iB^i* 



Set the index taking turns, that is, 71 "fn+2 and 7i + l — > 7i,(n + l) > i > 1). Again from the definition of di{E[^^ n], 
we can verify easily under the above index taking turns, 

■"+2 Q-iE.,A "+2 g-iB^.t 



g^-^)'"d.(i?[7,n+l]) ^ g(-^)'"d.(i;[7,n+l])- ^^^'^ 
Therefore we obtain the conclusion 




This implies that we finish the proof of expression (|103p . Substitute it into Eg. ([55)1 . we immediately obtain the same 
expression as our solution (|44)) . 

From the statement above, it can be seen that our solution actually finishes the task to solve the recurrence equation 
of bi^{t) by using our new method. Although the recurrence equation of 6!^(i) can be solved by integral in principle, 
it is indeed not easy if one does not know our identity (|30p and relevant relations. To our knowledge, the general 
term form has not ever been obtained up to now. This is an obvious difficulty to attempt absorbing high order 
approximations. However, our solution has done it, and so it will make some calculations of perturbation theory more 
convenient and more accurate. 

However, we pay for the price to require that H is not explicitly dependent on time. Obviously our solution is 
written as a symmetric and whole form. Moreover, our method gives up, at least in the sense of formalization, the 
requirement that V is small enough compared wit Hq, and our proof is more general and stricter. In addition, we 
do not need to consider the cases that the perturbing potential is switched at the initial and final time, but the 
perturbing potential is not explicitly dependent of time. Since we obtain the general term, our solution must have 
more applications, it is more efficient and more accurate for the practical applications, because we can selectively 
absorb the partial contributions from the high order even all order approximations. This can be called as the improved 
forms of perturbed solutions, which will be given in Ref. Q. 



VI. DISCUSSION AND CONCLUSION 



First of all, we would like to point out that our general and explicit solution (|44|) or its particular forms (|53|54p of 
the Schrodinger equation in a general quantum system independent of time is explicit and exact in form because all 
order approximations of perturbation not only are completely included but also are clearly expressed, although it is an 
infinite series. Our exact solution is in c-number function rather than operator form. This means that it can inherit 
the same advantage as the Feynman path integral expression. Moreover, our exact solution is a power series of the 
perturbing Hamiltonian like as the Dyson series in the interaction picture. This implies that the cut-off approximation 
of perturbation can be made based on the needed precision of the problems. It must be emphasized that when 
applying our exact solution to a concrete quantum system, all we need to do for our exact solution is the calculations 
of perturbing Hamiltonian matrix and the coefficient function limitations. Here, the "perturbing Hamiltonian matrix" 
refers to the representation matrix of the perturbing Hamiltonian in the unperturbed Hamiltonian representation. 
Thus, such calculations are definitely much easier than the calculations of path integral and Dyson series for the most 
systems. On the other hand, because the general term is explicitly given in our exact solution series, we can more 
conveniently consider the contributions from the high order even all order approximations in order to improve the 
precision than doing the same thing via the Dyson series. 

Although we have obtained the exact solution in general quantum systems, this does not mean that the perturbation 
theory is unnecessary because our exact solution is still an infinite power series of perturbation. Our solution is called 
"exact one" in the sense including the contributions from all order approximations of the perturbing Hamiltonian. In 
practise, if we do not intend to apply our solution to investigations of the formal theory of quantum mechanics, we 
need to cut off our exact solution series to some given order approximation in the calculations of concrete problems. 
Perhaps, one argues that our exact solution will back to the usual perturbation theory, and it is, at most, an explicit 
form that can bring out the efficiency amelioration. Nevertheless, the case is not so. Such a view, in fact, ignores the 
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significance of the general term in an infinite series, and forgets the technologies to deal with an infinite series in the 
present mathematics and physics. From our point of view, since the general term is known, we can systematically 
and reasonably absorb the partial contributions from some high order even all order approximations just like one has 
done in quantum field theory via summation over a series of different order but similar feature Feynman figures. In 
our paper [6], we will fixed this problem via proposing and developing an improved scheme of perturbation theory. 

In order to account for what is more in our solution, and reveal the relations and differences between our solution 
and the existed method, we compare our solution with the usual perturbation theory. We find their consistency. In 
fact, our solution has finished the task to calculate the expanding coefhcients of final state in Hq representation and 
obtain the general term up to any order using our own method. But the usual perturbation theory only carries out 
this task from some given order approximation to the next order approximation step by step. In a sentence, more 
explicit and general feature of our exact solution can lead to more physical conclusions and applications. In addition, 
by comparison with the existed perturbation theory, we actually verify the validity of our exact solution. 

It is worth pointing out that our solution, different from the time-dependent perturbation theory, presents the 
explicit solution of recurrence equation of expanding coefficients of final state in Hq representation, but we pay for 
the price that H is not explicitly dependent on the time. In short, there is gain and there is lose. 

It is worthy noticing that there are the apparent divergences in the expression of our exact solution, in special, 
when the degeneracy happens. These apparent divergences can be easily eliminated via the limitation calculations. 
The related discussions are arranged in our paper 6] . 

Therefore, we can say that our exact solution is more explicit than the usual non-perturbed solution of the 
Schrodinger equation because the general term is given in the c-number function form, and we can say our solu- 
tion is more general than the usual perturbation theory because our deducing methods give up some preconditions 
used in the usual scheme. However, our exact solution can not be used to the time-dependent system at present. 
Just as its explicit and general features, our exact solution not only has the mathematical delicateness, but also can 
contain more physical content, obtain the more efhciency and higher precision and result in new applications, which 
have been revealed in our serial studies 0, 01 . 

Our physical idea is a combination of Feynman path integral spirit and Dyson series kernel. Hence, we first find 
an general and explicit expression of the time evolution operator that not only is a c-number function but also a 
power series of perturbation including all order approximations. From our point of view, it is one of the most main 
results in our method. It must be emphasized that the study on the time evolution operator plays a central role in 
quantum dynamics and perturbation theory. Because the universal significance of our general and explicit expression 
of time evolution operator, we wish that it will have more applications in quantum theory. Besides the perturbation 
theory and open system dynamics, it is more interesting to apply our exact solution to the other formalization study 
of quantum dynamics in order to further and more powerfully show the advantages of our exact solution. 

In the process of deducing our exact solution, we obtain an expansion formula of operator binomials power. At our 
knowledge, this formula is first proposed and strictly proved. Besides its theoretical value in mathematics, we are 
sure it is interesting and important for expressing some useful operator formula in quantum physics. In addition, we 
prove an identity of fraction function. It should be interesting in mathematics. Perhaps, it has other applications to 
be expected finding. 

We can see that in our solution, every expanding coefhcient (amplitude) before the basis vector \^'') has some 
closed time evolution factors with exponential form, and includes all order approximations so that we can clearly 
understand the dynamical behavior of quantum systems. Although our solution is obtained in Hq representation, its 
form in other representations can be given by the representation transformation and the development factors with 
time t in the expanding coefficients (amplitude) do not change. 

Because we obtain the general term form of time evolution of quantum state, it provides the probability considering 
the partial contributions from the high order even all of order approximations. For example, we obtain the revised 
Fermi's Golden Rule [l3|. This can be seen in our recent work 

Our general and explicit solution can be successfully applied to o pen systems to obtain the master equation and 
exact solution. For example, we obtain the Redfield master equation [111 ] without using Born-Markov approximation. 
This can be seen in our recent work ,7} . 

From the features of our solution, we believe that it will have interesting applications in the calculation of en- 
tanglement dynamics and dccoherence process as well as the other physical quantities dependent on the expanding 
coefficients. Of course, our solution is an exact one, its advantages and features can not be fully revealed only via the 
perturbative method. 

In summary, our results can be thought of as theoretical developments of quantum dynamics, and are helpful 
for understanding the dynamical behavior and related subjects of general quantum systems in both theory and 
application. Together with our perturbation theory and open system dynamics Q they can finally form the 
foundation of theoretical formulism of quantum mechanics in general quantum systems. Further study on quantum 
mechanics of general quantum systems is on progressing. 
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APPENDIX A: THE PROOF OF OUR IDENTITY 



In this appendix, we would like to prove our identity ((30|) . For simplicity in notation and universality, we replace 
the variables E^^ by Xi as well as E['y^l] by x[l]. It is clear that the common denominator D{x[n]) in the above 
expression (fSO]) (the index I is replaced by n) reads 



while the i-th numerator is 



and the total numerator N[x[n],K) is 



D{x[n])=\{ 



ni{x\n\,K) 



Y\{x^-xj+i) 



D{x[n]) ^ 
di{x[n\y 



n+1 



N{x[n],K] ^^ni{x[n],K). 

i=l 

In order to simplify our notation, we denote ni{x[rL\) = ni{x[rL\, 0). Again, introducing a new vector 

2;fW = {a;2,a;3, • • ■ ,a;„,x„+i}. 



X^[n] = {Xi,- ■ ■ ,X^-l,X^J^l,■ ■ ■ ,Xn+l}, 

) 

+1 



•^n+il'^] {Xl^X2,^^^ , X ^ — I , X u) ' 



(Al) 



(A2) 



(A3) 



(A4) 
(A5) 
(A6) 



Obviously, xf{x[n]) with n components is obtained by deleting the i-th component from x[n]. From the definition of 
nj(x[n]), it follows that 

n,{x[n])^D{xf[n]). (A7) 

Without loss of generality, for an arbitrary given i, we always can rewrite — y[n — 1] = {yi, y2, • ■ ■ , y-n-i, Un} 

and consider the general expression of D{y[n — 1]). It is easy to verify that 



D{y[l])^-y2ni{y[l])+yin2{y[l]), 
D{y[2])^y2y3ni{y[2]) - yiy3n2{y[2]) + yiy2n3{y[2]). 

Thus, by mathematical induction, we first assume that for n > 2, 

n 

D{y[n - 1]) = ^(-l)(*-i)+("-iV,(y[n - l])n,(y[n - 1]), 

where we have defined 



i=l 



(A8) 
(A9) 



(AlO) 



(All) 



i=i 



As above, we have verified that the expression (jA10|) is valid for n = 2,3. Then, we need to prove that the following 
expression 



n+l 



D{y[n]) = J2i-lf-'^+yMn])n^iy[n]) 



(A12) 
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is correct. 

To our purpose, we start from the proof of a conclusion of the precondition (|A10[) as the foUowing: 



n+l 



^(-irv(j;W)-o. 



(A13) 



According to the relation (jA7p and substituting the precondition (jAlOp . we have 



ji+i 



n+l 



n+l 



1=1 



p,(yfN)n,(y,^[n]) 



n n+l 



Based on the definitions of pi{x[n]) and ni{x[n\), we find that 

p.(yf N) - 



P^-liyf[n]) (Ifz>j) 



n»-i(yf[n]) (Ifi>j) 



7i,(yf+Jn]) (Ifz<j) 



(A14) 

(A15) 
(A16) 



Thus, the right side of Eq. (|A14P becomes 



n n+l 

(-1)""' E E(-i)'^'?'^(2/^''N)n,(yf [n]) 
j=i i=i 

n n+l 

= (-i)""'E E i-^r^'pAy?[n])n,{y?W]) 

3 = 1 i=j+l 
n j 

j=i i=i 

n n+l 

= (-i)""'E E (-i)"+^K-i(2/fN)«.-i(yfM) 

n j 

+(-ir"'EE(-i)'^'p»(^f+iN)n.(yf+iN)- 



(A17) 



Setting that i — 1 i for the first term and j + 1 — » j for the second term in the right side of the above equation, we 
obtain 

n n+l 

j=i i=i 

n n 

= (-i)"-'EE(-i)^^'"'^'»(yfw)"^«(2^fN) 

+(-i)""'EE(-i)'^'"'^''(2^fN)«.(2/fM) 

j=2 ^=l 
n n 

= (-i)"EE(-i)^^'^''(yfW)"«(2^fw) 

n+l 

+ (-l)"5](-l)^+("+iV.(y^+iM)n,(y,f+iM). 



(A18) 
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Again, setting that i ^ j in the right side of the above equation gives 

n n+1 

j=i i=i 



Thus, it imphcs that 



n n+1 



(-ir"'EE(-i)^'''^'^(2^^''NK(yf[n]) = o. 

3 = 1 i=l 

From the relation (|A14[) it follows that the identity (jA13p is valid. 

Now, let us back to the proof of the expression (|A12[) . Since the definition of D{y[n]), we can rewrite it as 

D{y[n]) = D{y[n-l])d^+Mn]). 

Substituting our precondition (|A10p yields 

n 

D{y[n]) = 5^(-l)('-i)+("-i)p.(y[n-l]K(y[n-l])d„+i(yM). 

Note that 



(A19) 



(A20) 



i=l 



ni{y[n - l])dn+i{y[n]) = n^{y[n]){yi - yn+i), 
Pz{y[n - l])yn+i = Pi{y[n]), 
Pi{y[n-l])yi = pn+i{y[n]), 



(A21) 



(A22) 



(A23) 
(A24) 
(A25) 



we have that 



^(?iW)=E(-i)^'"'^^>'(2^N)n.(?iW) + E(-i)^'"'^^'""'Wi(2/N)n.(yN) 

n r n ' 

n 

= E(-i)^^"'^^>^(yW)"'(yW) + (-i)"-V+i(yN) [-(-i)"nn+i(yM)] . 

In the last equality we have used the conclusion (|A13p of our precondition, that is 

n 

Thus, Eq. (|A26|) becomes 

n+1 



(A26) 



(A27) 



(A28) 



The needed expression (jA12p is proved by mathematical induction. That is, we have proved that for any n > 1, the 
expression (|XT2l) is valid. 

Since our proof of the conclusion (jAlSp of our precondition is independent of ti (?! > 1), we can, in the same way, 
prove that the identity (|A13[) is correct for any n > 1. 

Now, let us prove our identity ([30)1 . Obviously, when X = we have 



n+1 



di{x[n]) D{x[n]) 



'n + 1 



= 0, 



(A29) 
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where we have used the fact the identity (jA13|) is vahd for any n > 1. Furthermore, we extend the definition domain 
of C^(x[n]) from K > n to K > 0, and stiU write its form as 

n+1 

Q-(xN) = g(-ir-^. (A30) 

Obviously, Eq. (jA29ll means 

CO(x[n])-0, (n>l). (A31) 
In order to consider the cases when if ^ 0, by using of di{x[n]){xi ~ Xn+2) = di{x[n + 1]) {i < n + 1), we obtain 



i— 1 

n+1 K / _ \ 

__ l)'~^ ^' I ^n+2 \ 



i=l 



1]) 



K+l n+1 X 

- CSHx|„ + ll,-,-l,-.-.5--±i-^-.„«g(-l)--'5^^ (A32) 

It follows the recurrence equation as the following 

C^{x[n]) = C^+,\x[n + 1]) - x„+2C,f+i(x[n + 1]). (A33) 

It implies that since C°(x[n]) ~ for any n > 1, then C\j^i{x[ii + 1]) = for any n > 1 or C\{x[n]) — for any 
n > 2; since C^(x[n]) = for any n > 2, then C'^_^i{x[n + 1]) = for any n > 2 or C'^{x[n]) = for any n > 3; ■ • • , 
since C^(x[n]) = for any n > (fc + 1) (fc > 0), then C'^Xii^l^ + 1]) = for any n > (fc + 1) or C^+^{x[n]) = for 
any n > (fc + 2); • • ■ . In fact, the mathematical induction tells us this result. Obviously 

C^{x[n]) = 0, (If < if < n). (A34) 

Taking if = n in Ea. (|A33p and using Eq. (|A34p . we have 

Cl{x[n])^ClXMn + l])^l, (n>l), (A35) 

where we have used the fact that Cj^(x[l]) = 1. Therefore, the proof of our identity ([50]) is finally finished. 
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